Topological insulator in a Bi-Bi 2 Se 3 infinitely adaptive superlattice phase 
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We report spin- and angle-resolved photoemission studies of a topological insulator from the 
infinitely adaptive series between elemental Bi and Bi 2 Se3. The compound, based on Bi4Se3, is a 1:1 
natural superlattice of alternating Bi 2 layers and Bi2Se3 layers; the inclusion of S allows the growth 
of large crystals, with the formula Bi4Se2.6S0.4- The crystals cleave along the interfaces between 
the Bi2 and Bi2Se3 layers, with the surfaces obtained having alternating Bi or Se termination. The 
resulting terraces, observed by photoemission electron microscopy, create avenues suitable for the 
study of one-dimensional topological physics. The electronic structure, determined by spin- and 
angle- resolved photoemission spectroscopy, shows the existence of a surface state that forms a 
large, hexagonally shaped Fermi surface around the F point of the surface Brillouin zone, with the 
spin structure indicating that this material is a topological insulator. 



PACS numbers: 74.25.Kc, 71.18.+y, 74.10,+v 

Three-dimensional (3D) topological insulators (TIs) 
have Dirac-like surface states in which the spin of the 
electron is locked perpendicular to its momentum in a 
chiral spin structure where electrons with opposite mo- 
menta have opposite spins [THE]. The number of exper- 
imentally identified 3D TIs has grown since the initial 
discovery, now including the Bii-^Sb^ alloy [3], Bi2Se3 
and Bi2Te3 and their mixed anion variants [6] \7\ I9HTT], 
GeBi 2 Te 4 and its PbSe-Bi 2 Se 3 variants [SHE], TlBiSe 2 
[131 [33], and SnTe [12]. Of these materials, Bi 2 Se 3 has 
been highlighted as a model system due to its large bulk 
gap and simple surface electronic structure - a single 
Dirac-like topological surface state (TSS) whose Dirac 
point resides inside the gap, essentially fully spin polar- 
ized and robust against the elastic scattering and scat- 
tering on phonons [4] [HJ [T5HT9] . 

Bi 2 Se3 is actually only one member of a family of com- 
pounds, called an "infinitely adaptive series", between 
pure Bi metal and Bi 2 Se3, made by stacking integer 
numbers of bismuth bi-layers and Bi 2 Se3 layers, in com- 
pounds of the type (Bi 2 ) m (Bi 2 Se3) n [20]. Given the topo- 
logical properties of the surface of Bi 2 Se3, and the fact 
that the bismuth bi-layer has been predicted to be a two- 
dimensional quantum spin Hall system [21 , it is of inter- 
est to determine whether the superlattice material also 
has topological surface states. Thus far, however, these 
materials have been studied only structurally as small 
crystals, or in polycrystalline form [20]. Here we report 
the successful growth of large single crystals of the sim- 



plest member of the infinitely adaptive series, i.e. Bi4Se3 
(m=n=l), stabilized for crystal growth by partial S sub- 
stitution for Se, and present spin- and angle-resolved pho- 
toemission spectroscopy (SARPES) studies of its surface 
electronic structure. Our measurements reveal a hexago- 
nally shaped surface state whose spin structure indicates 
that Bi4Se3 is a new topological insulator. The com- 
pound is the natural bulk analog of recently described 
artificial material made by Bi adsorption on the surface of 
Bi 2 Te3 [22] . Photoemission electron microscopy (PEEM) 
reveals that the material cleaves with mixed Bi-layer and 
Se-layer terminations of the exposed terraces. This pro- 
hibits the association of the TSS with either the Bi or 
Se terraces by (S)ARPES, but it shows that this mate- 
rial offers the unique opportunity for studying novel ID 
topological physics near the terrace steps by nanometer- 
scale probes. 

High purity elemental Bi (99.999%), Se (99.999%), 
and S (99.999%) were used for the Bi4Se 2 .6So.4 crystal 
growth, with a starting formula of Bi4.5Se 2 S; the excess 
bismuth was used as a flux and the S greatly enhances 
the selective crystal growth of the 4:3 phase. The pre- 
cleaned elements were sealed in clean quartz ampoules, 
heated to 850 °C for 1 day, and cooling over the pe- 
riod of 1 hour to 700 °C. The crystal growth involved 
cooling from 700 °C to 550 °C over a period of 25 h, 
followed by centrifuging. Single crystals up to 1 centime- 
ter size were obtained, with high-luster cleavage faces. 
Single crystal X-ray diffraction data was collected on a 
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FIG. 1: (Color on line) Crystal structure (a) and Photoemis- 
sion Electron Microscopy (PEEM) images of a Bi4Se2.6So.4 
surface, taken at the Bi 5d (b) and Se 3d (c) core levels with 
100 eV photons. Field of view in the PEEM images is 20x20 
/im. Bright (dark) contours reflect regions of high (low) pho- 
toemission intensities. 



Bruker APEX II diffractometer using Mo Ka radiation 
at room temperature. Unit cell determination and refine- 
ment, and data integration were performed with Bruker 
APEX2 software. The crystal structure was determined 
using SHELXL-97 implemented through WinGX. 

The ARPES experiments were carried out on a Sci- 
enta SES-100 electron spectrometer at beamline 12.0.1 of 
the Advanced Light Source. The spectra were recorded 
at photon energies ranging from 35 to 100 eV, as indi- 
cated in the text, with a combined instrumental energy 
resolution of ~ 12 meV and an angular resolution bet- 
ter than ±0.07°. Samples were cleaved at low temper- 
ature (15-20 K) under ultra-high vacuum (UHV) con- 
ditions. The temperature was measured using a silicon 
sensor mounted near the sample. The S ARPES experi- 
ments were performed at beamline U5UA at the National 
Synchrotron Light Source using a Scienta R4000 electron 
analyzer coupled to a pair of mini-Mott spin polarime- 
ters that enable measurements of both the in-plane and 
out-of plane spin polarizations. The spin-resolved data 
were recorded at 58 eV photon energy, with the sam- 
ple kept at ~ 300 K. Energy and angle resolution were 
approximately 40 meV and 0.5°, respectively. Photoe- 
mission Electron Microscopy (PEEM) experiments were 
carried out at the XPEEM/LEEM end-station at the Na- 
tional Synchrotron Light Source beamline U5UA at room 
temperature, with 100 eV photons using an Elmitec SPE- 
LEEM III microscope. All the samples were cut from the 
same bulk piece and cleaved and measured in ultrahigh 
vacuum conditions (base pressure better than 2 x 10 -9 Pa 
in the ARPES and SARPES chambers and better than 
2 x 10~ 8 Pa in the PEEM chamber). 

Electronic structure calculations where performed in 
the framework of density functional theory (DFT) us- 



TABLE I: Structural Characterization of Bi4Se2.56S0.44. 
Space group R3(No. 148), a = 4.2466(8) A, c = 39.706(7) 
A, 329 unique reflections, Ri (all reflections) = 0.0457, WR2 
= 0.0938, G00F = 1.109. 
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ing the wien2k [23] code with a full-potential linearized 
augmented plane-wave and local orbitals [FP-LAPW + 
lo] basis [24-26 together with the PBE parameterization 
[27] of the GGA for the exchange-correlation functional. 
The plane wave cut-off parameter Rmt^max was set 
to 7 and the irreducible Brillouin zone (BZ) was sampled 
by 288 k-points. Spin orbit coupling (SOC) was included. 
For Bi4Se3 the reported crystal structure [28] was used. 
To estimate the effect of the partial S substitution for 
Se, electronic structure calculations were also performed 
for Bi 4 Se2S. For this material a fully ordered structure 
was taken (with S fully occupying the middle layer in the 
quintuple layer part of the structure) and the lattice con- 
stants were interpolated linearly, according to Vegards's 
law [29], using our data for the structure of Bi4Se2.6S0.4- 

Figure [I] shows a schematic view of the crystal struc- 
ture of Bi 4 Se2.6So.4 and PEEM images of a cleaved sur- 
face. The details of the structural characterization are 
given in Table 1 ; the structure determination was used to 
precisely define the crystal composition. The two PEEM 
micrographs show the same spot on the surface recorded 
at two different kinetic energies of photo-electrons, corre- 
sponding to the Bi bd and Se 3d core levels. The inversion 
of contrast between Bi 5d and Se 3d micrographs indi- 
cates that the surface consists of alternating Bi2 and Se 
terminated terraces. The terrace width is of the order of 
a few /im. As the photon spot sizes in our ARPES and 
SARPES experiments are around 50 /im and 0.5 mm, re- 
spectively, the photoelectron characterization represents 
an average over the two types of terraces. 

In Figure [2] we compare the experimentally determined 
electronic structures of Bi4Se2.6So.4 and Bi2Se3. The 
Bismuth core level (5d) spectra (Fig. [2ja)) clearly indi- 
cate two different chemical environments for Bi atoms in 
Bi 4 Se2.6So.4' one is similar to the environment in Bi 2 Se3, 
whereas the other indicates the more metallic character of 
the Bi2 layers. The wide-energy range valence band spec- 
tra, shown in Figure[2](b) and (c) are generally similar for 
the two materials, with a nearly rigid-band down-shift of 
~ 0.4 eV of the Bi 4 Se2.6So.4 spectrum relative to Bi 2 Se3 
yielding a good general correspondence of the bands. The 
high intensity feature in the Bi4Se2.6So.4 spectrum that 
crosses the Fermi level is the central point of this study 
- this is the TSS, similar to the one observed in Bi2Se3, 
but also down-shifted, with the Dirac point around 0.54 
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FIG. 2: (Color on line) Photoemission Spectra from 
Bi4Se2.6So.4 and Bi2Se3. (a) Bi 5d core-levels in Bi4Se2.6So.4 
(blue) and Bi 2 Se3 (red). Double peak structure of each spin- 
orbit split state indicates two chemically different bismuth 
sites in the Bi4Se2.6So.4 crystal. Inset shows the SBZ with the 
Bi4Se2.6So.4 Fermi surface measured with 59 eV photons, (b) 
Valence band along the Y — M line in the SBZ for Bi4Se2.6So.4 
and (c) for Bi2Se3. 



eV below the Fermi level. This down-shift of both the va- 
lence bands and the TSS suggests that the major effect 
of the Bi2 layers in this natural superlattice phase is to 
donate electrons to the Bi2Se3 quintuple layers. A wide 
/c-space survey indicates that the TSS forms a hexago- 
nally shaped Fermi surface at the center of the surface 
Brillouin zone (SBZ) (Fig. ga), inset). A new feature 
in the valence band spectrum in Fig. [2|b) is the "hole- 
like" state that crosses with the TSS at —0.15 eV, 
k x « iO.lSA- 1 . This is a bulk state, as will be discussed 
further with reference to Fig. [3j which does not have an 
obvious counterpart in the Bi2Se3 spectrum. 

In Figure [3j we show the detailed spin- and k- resolved 
electronic structure of Bi4Se2.6S0.4- Panels (a) and (b) 
show the photon energy dependence of the ARPES in- 
tensity at the Fermi level along the V — K line and at 
the center of the SBZ, respectively, for photon energies 
ranging from 50 to 84 eV in 2 eV steps. The intensity 
maps are converted into /c-space by using the inner po- 
tential V = 10 eV in the approximation for the out-of- 
plane momentum, k z = j-^/2m(EkCos 2 (0) + V), where 
Ek is the kinetic energy of a photoelectron. The straight 
features at k y = ±0.173 A -1 in Fig. j3^a) represent the 
Fermi surface of the TSS. The bulk hole state previously 
mentioned is visible here for negative k y as a slightly dis- 
persing (with k z ) feature, intertwined around the TSS. 
The straight vertical line of high intensity in (b) corre- 
sponds to the Dirac point of the TSS. Slightly lower in 
energy, another bulk state can be seen as dispersing with 
k z , with the minimum around k z = 4.35 A -1 . Both bulk 
features have the periodicity of « 0.5 A -1 , suggesting 
that the relevant periodicity in real space in the z direc- 
tion is « 12.5 A. This is close to the thickness of one of 
the Bi2-Bi2Se3 composite layers that repeats three times 
to form the rhombohedral unit cell (Fig. [IJa)). In panels 



(c) and (f), we also show the in-plane ARPES intensity 
at the Fermi level for two photon energies, hv = 58 and 
72 eV, corresponding to k z = 4.06 and 4.43 A , re- 
spectively. The hexagonally shaped contour represents 
the Fermi surface of the TSS. The dispersion of the sur- 
face state measured at hv = 58 eV along two high sym- 
metry directions in the SBZ is shown in panels (d) and 
(e). At hv = 72 eV, the intensity at the Fermi level is 
strongly modulated by the nearly overlapping bulk state 
with the three- fold symmetry (Fig. [3^f ) ) . The area en- 
closed by the TSS' Fermi surface is almost an order of 
magnitude larger than that in pristine Bi 2 Se3, reflecting 
a surface electron concentration of ~ 8.7 x 10 12 cm -2 . 
Such high surface electron doping can also be obtained 
on Bi2Se3 through adsorption of various metals on it's 
surface, producing a very similar hexagonal warping of 
the TSS' Fermi surface [18]. 

Figure j3^g) shows the spin-resolved spectra with po- 
larization defined in the y direction for points A and 
B from Fig. [3^c). The peak at the Fermi level, corre- 
sponding to the TSS, shows a significant spin polariza- 
tion (~ 30%) that reverses sign when the momentum is 
reversed (k x —> —k x ), an indication of the chiral spin 
structure. It is interesting to note that the chirality is 
opposite to that observed in Bi 2 Se3 |8]. The peak at 
around —0.5 eV also shows a polarization, but without 
the sign-reversal, an effect also observed in the bulk va- 
lence band in Bi2Se3 [15]. The polarization of the TSS is 
much smaller in Bi4Se2.6So.4 than in Bi2Se3 [8 , probably 
due to the partial overlap with the bulk states. 

Figure [4] shows a comparison of the calculated bulk 
band structures for Bi2Se3 and Bi4Se3; SOC has a signif- 
icant effect on the band structure of Bi 4 Se3. The lower 
red colored band in figure [4] (b), for example, has shifted 
by about 1 eV on the inclusion of SOC. The calculations 
predict Bi4Se3 to be a semi metal with hole pockets at 
T and Z and an electron pocket at F. The hole pockets 
at T and Z project to the f point in the SBZ and are 
seen in the ARPES measurement. The electron pocket 
predicted at the bulk F point is not observed by ARPES, 
suggesting that it lies higher in energy than is predicted 
by the calculations. Comparison of the band structures 
allows us to identify bands closely above and below in 
Bi 4 Se3 as being Bi 2 Se3-like bands that have been modi- 
fied through interactions with the Bi bilayers; these are 
shown as red in the figure. Bands appearing directly 
around Ep (blue) in Bi4Se3 are primarily Bi2 derived 
states modified by hybridization with the Bi2Se3 layers. 

There are two continuous gaps in the electronic struc- 
ture for Bi4Se3, labeled as Gap 1 and Gap 2 in Figure [4] 
Counting all bands to Gap 1 gives, for the time-reversal- 
invariant momenta (TRIM) points of the bulk BZ [30] , 
positive parity for T, Z, and L, and negative parity for 
F. Therefore in Bi4Se3 we expect a TSS cone at the pro- 
jection of the F point in the SBZ. According to the band 
structure calculation, the cone would be about 0.5 eV 
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FIG. 3: (Color on line) Spin-resolved electronic structure of Bi4Se2.6S0.4- (a) ARPES intensity at the Fermi level as a 
function of the in-plane (k y along the Y — K line) and out-of plane (k z ) momentum. The two vertical lines of high intensity 
at k y = zLkF ~ ±0.173 A correspond to the Fermi surface of the surface state from Fig. 2(b). (b) Intensity at the Y point 



The vertical line of high intensity at E= 0.54 eV corresponds to the Dirac 



(fee = k y = 0) of the SBZ as a function of k z 

point of the surface state, (c) In-plane Fermi surface measured at k z = 4.06 A x and (f) at k z — 4.43 A x , corresponding to 
the red and blue lines in panel (a), respectively. The circles A and B in (c) indicate the /c-positions where the spin-resolved 
measurements were performed, while arrows represent the in-plane spin polarization. Dispersion of the surface state along the 
r — M (d) and Y — K (e) momentum lines in the SBZ, measured at k z — 4.06 A . (g) S ARPES spectra taken at A and B 
points in the SBZ, as indicated in (c). Red (blue) spectra represent spin-up (spin-down) components in the y direction. TSS 
is indicated by arrow. 



below the Fermi energy (for Bi4Se2S we expect a second 
cone at the projection of the Z point, which also has nega- 
tive parity, but since this would be above the Fermi level 
it would not be seen by ARPES). Counting the bands 
until Gap 2 yields, for both Bi 4 Se 3 and Bi 4 Se 2 S, TRIM 
points Z, M, and L with positive party and Y with nega- 
tive parity. This would yield a TSS cone at 0.5 eV below 
the Fermi level at Y. Note that counting to Gap 2 is 
equivalent to counting to the lower red band belonging 
to Bi2Se3; we observe the same parity at each TRIM as is 
seen in pure Bi2Se3. When counting to Gap 1 one counts 
to a band belonging to the Bi bi-layer. The second cone 
at F might belong to the Bi 2 , which was predicted to be 
a 2D quantum spin Hall system [21]. The crystal cleaves 
to expose Bi 2 and Bi2Se3 Se-terminated terraces, neither 
of which are at a bulk inversion center. In this case the 
projection of the bulk BZ points to the surface BZ is less 
certain, but since the only TRIM at the surface BZ are at 
f and M, the TSS must be at either of those two points. 
In ARPES we see one TSS cone at f and no cone at M 
(Fig. [2|a), inset). Where these surface states reside, i.e. 
on the Bi 2 terraces, the Bi2Se3 terraces, or both, and 
their relation to the bulk band structure, should be the 
subject of further study on this material. 

Bi4Se2.6So.4 opens several interesting possibilities for 
studying the role of dimensionality in topological sys- 
tems: as Fig. [3] shows, there is a 2D n-type TSS, nearly 
overlapping with a p-type bulk (3D) Fermi surface. In 



addition, PEEM images from Fig {T] show alternating Bi2 
- Se terminations, each kind of surface no doubt display- 
ing a different surface electronic structure. Each step 
would then represent a boundary (or a junction) between 
two 2D (topological-'trivial' or topological-topological, n- 
n or p-n) electronic "gasses". This offers the possibility 
of studying topological phenomena ranging from ID to 
2D to 3D in a single system. The current ARPES exper- 
iments lack the necessary spatial resolution to observe 
such effects, but the scanning tunneling microscopy, for 
example, could provide information about these phenom- 
ena. A recent study of deposited Bi bi-layer on Bi2Te3 
[22] has indicated that similar phenomena might be at 
play in that system: with depositing bismuth, a hole-like 
state is added to the Bi2Te3 bands. Although the au- 
thors consider the Bi2 layer to be homogeneous, it might 
be that terraces with mixed termination are formed in 
that case as well, creating a situation that is similar to 
what we observe in Bi4Se2.6S0.4- 
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FIG. 4: (Color on line) Comparison of the bulk electronic 
band structures for (a) Bi2Se3 and (b) Bi4Se3. Colored bands 
are described in the text. 
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